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Abstract 

Parton distributions of pseudoscalar tt,K and ij mesons obtained within the NJL 
model using the Pauli-Villars regularization method are analyzed in terms of LO 
and NLO evolution, and the valence sea quark and gluon parton distributions for 
the pion are obtained at = 4GeV^ and compared to existing parametrizations 
at that scale. Surprisingly, the NLO order effects turn out to be small compared to 
the LO ones. The valence distributions are in good agreement with experimental 
analyses, but the gluon and sea distributions come out to be softer in the high-x 
region and harder in the low-x region than the experimental analyses suggest. 
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1 Introduction 

The study of structure functions of hadrons in the Bjorken limit and high enough Q^, 
as probed in inclusive deep inelastic scattering (DIS), is traditionally considered the do- 
main of perturbative QCD since the running coupling constant becomes small [|^. 
Present day QCD leading order (LO) and next to leading order (NLO) phenomenological 
calculations can relate leading twist contributions to structure functions among differ- 
ent momentum scales through the well known linear integro-differential Gribov-Lipatov- 
Altarelli-Parisi (GLAP) equations 0. This makes sense if is high enough so that only 
leading twist logarithmic corrections contribute and higher twist power-like corrections 
are negligible. To start with, some theoretical or experimental nonperturbative profile 
function is needed as initial condition for the GLAP equations. In the nucleon case, QCD 
scaling violations have been confirmed by relating experimental partonic distributions 
at many values, and many phenomenological parametrizations have been proposed 
[0, 1^, 1^, ^. Naturally, these parametrizations are under continuous update to incorpo- 
rate increasing information obtained from current experiments. The net result is that 
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uncertainties in the parton distribution functions of the nucleon include not only a large 
body of experimental data, but also theoretical NNLO or higher twist error estimates 
which provide a, perhaps inaccurate, but undoubtedly systematic description within a 
large region in the x, plane (see, e.g., the talks in Ref. 0). 

In contrast with the nucleon case, our present knowledge of parton distribution func- 
tions for other hadrons is rather poor. As suggested long ago it is possible to estimate 
distribution functions using constituent quark models to evaluate the low energy initial 
condition under the assumption that the gluon and sea content of hadrons vanish at 
the corresponding low energy resolution scale, and dynamically generate them by QCD 
evolution to higher scales. These estimates can then be used to test the sensitivity 
of various experiments to the distributions of interest. Recently, this approach has also 
been applied to generalized parton distributions (GPD's) p|, 0, |ll], 0, |13|, [1^, which are 
generalizations of the usual parton distributions considered in this work and are related 
via a sum rule to the elastic form factors. Unlike the usual parton distributions, the 
GPD's are not directly measurable as observables are always expressed in terms of them 
via a convolution formula. This has lead some [|15] to conclude that, at present, model 
calculations of the GPD's at a low scale are needed to assess the sensitivity of various 
observables to the GPD's. In practice, this requires an evolution of the GPD's calculated 
at the low energy scale to the scale relevant to the experiment. As a is rather large at 
the low energy scale, one must worry about the use of perturbative evolution to connect 
the low energy model with the high energy data. Thus, it seems prudent to test this 
procedure in a situation where data are available to compare with theory, i.e., the usual 
parton distributions. To test the validity of this approach, it is necessary to compare the 
LO and NLO results not only for the valence distributions, but also for the sea and gluon 
distributions. 

From a theoretical viewpoint, pseudoscalar mesons and specifically vr and K mesons 
are particularly distinguished hadrons since most of their low energy properties follow 
the patterns dictated by chiral symmetry. Actually, we do not expect to understand the 
properties of any hadron better than the pion, as Chiral Perturbation Theory suggests. 
By extension, one might think that the parton structure of a pion is the simplest one 
to consider provided chiral symmetry constraints, i.e., spontaneous and explicit chiral 
symmetry breaking, are properly incorporated. The recent work [0 clarifies this point 
regarding explicit chiral symmetry breaking; ChPT allows one to systematically compute 
chiral corrections to the moments of structure functions, but says nothing about the 
soft pion limit. Each moment corresponds to a undetermined low energy parameter 
which renormalizes a local operator. On the other hand, improved QCD sum rules have 
also been employed |]T7[ to determine the quark distribution functions of the pion in the 
intermediate x region, 0.15 < X < 0.7 at Q2 = 2GeV^ where the model is applicable. The 
absolute normalization becomes a problem since some ansatz for the distributions must 
be made outside this x— range. 

Among the quark models where spontaneous breaking of chiral symmetry plays a 
dominant role, the Nambu-Jona-Lasinio (NJL) model provides a particular example of a 
chiral quark model where a unified picture of vacuum, mesons and nucleons is achieved 



18]. Pseudoscalar mesons appear as quark-antiquark excitations of the spontaneously 
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broken vacuum. Several calculations of the pion structure functions within chiral quark 
loop models have been made and many different results for the initial conditions have 
been obtained. One important and tricky reason for the discrepancies lies in the use 
of different regularization procedures. As the bosonized version of the NJL model is 
similar to other quark-loop models of the pion (the vrgg coupling is 75-like), we think 
it of interest to briefly review them and comment on the main differences. The use 
of different regularizations might be regarded as an objection to the NJL model itself. 
However, not every regularization scheme can be considered acceptable. Actually, some 
of the quark-loop calculations violate some necessary conditions on the regularization. 
We argue in the following that in some cases one should blame the regularization scheme 
instead of the model. At a formal level, the process of going from the hadronic to the 
distribution function can be done by using the so-called quark-target scattering formula 
[pj[] . A large body of quark loop model calculations have been done making use of these 
ideas 
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28l |29|, pOl pll]. The problem of proceeding in that 



way is that the distribution function may turn out to be non-normalizable. 

Unfortunately, no first principle calculations of structure functions for pseudoscalar 
mesons are yet available with the exception of some standard lattice calculations of the 
lowest moments 0, albeit in the quenched approximation and subjected 

to well-known problems with chiral extrapolations. In addition, the reconstruction of 
the X dependent structure function via some inverse moments method is strongly biased 
in the intermediate and low x regions. The transverse lattice approach employed in 
Refs. I^G], ^ offers the possibility of directly computing structure functions in x-space. In 
any case, as one might expect from the quenched approximation, the lattice results provide 
a larger momentum fraction of valence quarks than those suggested by phenomenological 
|31 



analyses 



Although not as well determined as the nucleon, the parton structure of the pion 
has been analyzed on a phenomenological level ^8] and a simple parametrization at 
Q2 = 4GeV2 has been given. The valence quark distributions extracted in this work pHI 



from Drell-Yan experiments seem well determined, whereas the gluon distribution as 
obtained from the analysis of prompt-photon emission data is less well determined. 



On a phenomenological level, the constituent model proposed in Ref. ||4^ for the valence 
distributions of the pion has been further extended to the sea and gluon distributions |^ 
and the K/n valence up-quark ratio. In these calculations and in the recent update 
the required total valence momentum fraction in the pion at = 4GeV^ is taken to be 
the same as for the nucleon, (xV^) = (xVn) = 0.40, a bit below the value (x^-) = 0.47 of 
Ref. PS|. As different data sets have been fitted and different nucleon parton distributions 
have been used in the different analyses, it is not clear what to make of the differences. In 



addition, although Ref. |^ includes error estimates, the model analysis of Refs. ^| 
does not include them, and therefore it is not possible to know if the differences are 
significant. Let us note that the E615 experiment |^ suggests the valence density of 
the pion may be enhanced by about 20% compared to the proton, and a recent analysis 
4^ of the ZEUS di-jet data seem to favor the gluon distributions of Ref. ||3^. Thus, in 



determining the low-energy scale of our model, we use the valence momentum fraction 
found in Ref. p8[. Finally, we also compare with the K^/tt^ structure function ratio at 
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Q2 = 20GeV^ which has been measured using the Drell-Yan process 
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Remarks on Pion parton distribution functions in 
quark loop models 



In a previous work we found the structure function of the pion to be a constant 
function of x in the NJL model in the chiral hmit and in the leading order of a large Nc 
expansion. To get this result the use of a suitable regularization method was needed. A 



thorough study of regularization methods in the NJL model may be found in Ref. 
and we refer to that work for a more detailed description. By suitable we mean several 
desirable properties that should be incorporated, namely: 



• The connection between the forward Compton amplitude and the quark-target scat- 
tering amplitude is valid only for gauge invariant, finite amplitudes. For this reason, 
some gauge invariant regularization must be imposed on the Compton amplitude. 
Naive sharp cut-offs are not acceptable from this viewpoint. In addition, this way 
of proceeding represents a further advantage, since in the NJL model it is only 
known how to regularize closed quark loops. The quark target scattering amplitude 
corresponds to an open quark line. 

• The regularization must produce exact scaling in the Bjorken limit. The main reason 
is that this is the only way we know how to extract the leading, and eventually 
higher, twist contributions for which QCD evolution is known. This eliminates 
proper-time regularization, since it produces unrealistic scaling violations. 

• The regularization must also be able to work away from the chiral limit, but without 
spoiling the QCD anomaly. The former condition precludes a single Pauli-Villars 
subtraction. 

• The regularization should allow calculations in both Minkowski and Euclidean space, 
i.e., dispersion relations must be fulfilled. This turns out to be very convenient for 
DIS calculations, since cutting rules may be used. 

• The resulting distributions should satisfy the normalization condition and the mo- 
mentum sum rule. 



We found in Ref. 22 



that the Pauli-Villars with two subtractions fulfills the desired 
requirements. In addition, the Pauli-Villars scheme does not spoil the good description of 
other low energy hadronic properties found in the NJL model |1T8| , ^], fulfills dispersion 
relations [|^, and allows one to regularize the Dirac sea of the chiral soliton away from 

for definiteness, one gets in the chiral limit 



the chiral limit 26 



Taking vr^ 



d^+{x,Ql) = 9{x)9{l-x) 



(1) 
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The results for m,r 7^ are displayed for completeness in the Appendix. By construction, 
Eq.(|I|) is consistent with chiral symmetry. The result was obtained by several means 
within the NJL model either using Pauli-Villars regularization p2|, p6| on the virtual 



Compton amplitude or imposing a transverse cut-off p5| upon the quark-target ampli- 
tude. This result has been recently re-derived in a chiral quark model solving chiral 
Ward identities by using the so-called gauge technique |Q . The easiest way to understand 
Eq. (|^) is perhaps in terms of phase space arguments and point couphngs (i.e., constant 
matrix elements) For a massless pion this is justified since intermediate states in the 
quark-target amplitude have = mT^{l — x) and the low momentum components 
of Trqq matrix element dominate. Let us mention that Eq.(|l|) disagrees with other NJL 
calculations, due to the use of different regularizations. If the virtual Compton amplitude 
is used with a four-dimensional cut-off pOf or the quark-target amplitude is used with 



Lepage-Brodsky regularization [25], different shapes for the quark distributions are ob- 
tained. The null-plane |^T|] NJL model with sharp cut-off light-cone (LC) quantized 
NJL model and spectator model p4| calculations also produce different results. In all 
cases, the use of momentum dependent form factors or non-gauge invariant regularizations 
make the connection between Compton amplitude and quark-target amplitude doubtful 
and, furthermore, spoil normalization. The results based on a quark loop with momentum 



dependent quark masses ^ |3^ seem to produce a non-constant distribution. Recent 
calculations on the transverse lattice reveal either an almost flat structure very much 
resembling Eq. (|1|) at a scale = IGeV^ or a more bumped form 1^^. The reason for 
the discrepancy between these two transverse-lattice calculations is not obvious to us. 

In this paper we study within LO and NLO the parton content of pseudoscalar mesons, 
namely vr, K and rj including valence, gluon and sea distributions, thus extending our pre- 
vious work where only the initial conditions were presented and the LO evolution for 
the valence distributions. There, we analyzed the LO valence contribution and impressive 
agreement with SMRS PB| parametrization at = 4GeV^ was obtained. Encouraged 
by this success we extend our analysis to the sea and gluon distributions both in LO and 
NLO evolution. 



3 Numerical Results 



3.1 Momentum fraction analysis 



To perform the evolution, one must determine the scale Qo of the model. We determine 
this scale by fitting the valence quark momentum fraction at 4 GeV^. For definiteness, 
we take the running strong coupling constant at the Z mass, Mz = 91.12GeV, to be 
0.116 and evolve it down by exactly solving the differential equation 



a(M|) 



da 
Itt 



(3{a) 



-a 



a y 
4^/ 



+ 



(2) 



where t = \n{Q'^/Ql) and a = g^/{ATT). We take the number of active flavors to diminish 



by one unit each time a quark threshold is crossed, i.e., Nf{Q^) = J2u 
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Figure 1: Valence, gluon and sea distributions in the pion, vr'^, at = 4GeV^ in the NJL 



model compared with phenomenological analysis for the pion SMRS92 |^ and GRS99 
We take the valence momentum fraction {xV)-,^ = 0.47 at = 4GeV^. 



with TTib = 4.5GeV and rric = 2.0GeV. For Np = 3 one has /3o = 9, /3i = 64. This yields 
the value a(4GeV^) = 0.284. Below that scale we fix the number of flavors equal to three, 
since we consider evolution below charm threshold. The previous formula, Eq. (||), is used 
to transform the variable t into the variable a, by exactly^ solving the differential equation. 
Since we numerically perform the NLO evolution of the sea and gluon distributions, it is 
convenient to specify the initial at ti and numerically integrate Eq. (Q) to t/ to obtain 
aj. We note, however, that it is also possible to find an implicit solution for a/ in terms 
of ai, U and tj. Specifically, we find 

aj AttPo 47r/5o \l + ^a, ] 47r^ ^ " a, 

/3 //n2\ 1 

(3) 

Although we are not aware of an analytic solution for in terms of the other parameters, 
this equation may be solved numerically very quickly and accurately using Newton's 




^By 'exact' we mean solving the differential equation with a numerical accuracy much greater than 
the experimental uncertainty on a. 
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Figure 2: Valence, gluon and sea distributions in the kaon, K~^, at = 4GeV^ in the 
NJL modeL We take the total valence momentum fraction {xV)-,^ = 0.47 at = 4GeV^. 

method. For example, taking Qi = rric and = 0.284, and using this as the initial 
seed for aj, one obtains at least eight-significant digit accuracy for af after at most ten 
iterations all the way down to Q/ of 0.4 GeV. This form also enables one to determine at 
what scale aj diverges. For af ^ oo, we obtain 

Evidently, a/ diverges when the right-hand-side of the above equation vanishes, which 
happens at a scale of Qf ^ 0.365 GeV. 

The non-singlet momentum fraction satisfies the differential equation 

/3(a)^^^ = 72,7V5(«)V^2(«) (5) 
da 

where V2{a) can be any non-singlet quark distribution. Up to two loops one obtains the 
expansion p[ 



7.Ns(«)=7Ssf£)+7as(£) +••• (6) 
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Figure 3: up and strange quark and gluon LO and NLO distribution functions in the 
//-meson, at = 4GeV2 in the NJL model. 



To proceed further, we use the results from Ref. where it was found that at = 
4GeV^ valence quarks carry 47% of the total momentum fraction in the pion, e.g., for vr"*", 

{x (u^ -u^ + d^-d^^) = 0.47 at Q2 ^ ^Q^y2 _ ^-^-j 

Evolving downwards, we get that for ao = 1.89(1.487) valence quarks carry 100% of the 
total momentum in the pion in LO (NLO). 



3.2 Pion structure functions 



Having determined Qq of our model, we evolve the structure functions to = 4 GeV^ 



using the scheme presented in Ref. |£2[, which requires an analytical formula for the 
moments of the distribution function. In the chiral limit, = 0, the moments may 
be trivially computed. Away from it, 171^^ ^ , they can be expressed in terms of hy- 
pergeometric functions 2-^1, but it is more convenient, and just as accurate, to make a 
polynomial approximation in a;— space and then compute analytically the moments. For 
completeness, the result of such a fit for both u{x) and d{x) is presented in the appendix. 
For the 77, it is more convenient to make an expansion in terms of x(l — x), as is discussed 
in the appendix. In this work, we take mT^= 139.6 MeV, = 494 MeV, /^r = 93.3 MeV 
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Figure 4: Chiral corrections to the u-quark valence LO an NLO distribution functions 
at = 4GeV^ compared with phenomenological analysis for the pion SMRS92 |^8|] and 
GRS99 m. As suggested in Ref. m we take (xV)^ = 0.47 at = 4GeV^ 



and Mu = Md = 280 MeV, resulting in Ms= 527 MeV, A = 870 MeV and = 501 MeV 
(exp. 549 MeV). 

Our LO and NLO valence, sea and gluon distribution functions evolved from the 
quark model point, Ql, where the valence quarks carry all the momentum, to the point 

= 4GeV^ where gluon and sea distributions are dynamically generated, are shown in 
Fig. (0). They are compared to the phenomenological analysis of Refs. |^ and The 
remaining distributions trivially fulfill 

J^+(a;, Q^) = M^+(a;, Q^) u-^+ix^Q"^) = d^+ix^Q"^) Q^) = Q^) (8) 

as a consequence of our initial condition and properties of evolution. The LO valence 



result was already presented in our previous work ||2^. We see here that NLO evolution 
does not make a big difference, providing some confidence in perturbative evolution, even 
though the quark model point corresponds to a's larger than unity. Actually, it has been 
suggested that the natural expansion parameter for DIS is a/n, which in our case is about 
a half, a{Qo)/7r ~ 0.5. 

As can be seen in Fig. (^, the finite pion mass effects turn out to be rather small 
because chiral corrections to the initial condition are small within the model at the one 
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Figure 5: Dependence of the u-quark valence LO an NLO distribution functions at = 
4GeV^ on the momentum fraction at that scale, compared with phenomenological analysis 



for the pion SMRS92 |8[ which takes {xV)^ = 0.47 and GRS99 |g where {xV)^ = 0.40 
is used. 

loop level. While it is conceivable that pion loop effects could provide, as is frequently 
the case, some logarithmic enhancement to chiral corrections, it is a feature of GLAP 
evolution equations that upward evolution tends to wash out the differences in the initial 
condition. 

We finish our discussion on the pion parton distribution by comparing in Fig. (j^) the 
results obtained by taking either {xV)-,^ = 0.47 as suggested by the SMRS92 analysis 



3^ or {xV)t, = 0.40 as implied by the GRS99 parametrization [^. The sea and gluon 
distributions are not shown because their dependence on the momentum fraction is rather 
small. As can be deduced from the figure, the shape of the valence distribution is much 



better described if, as determined in Ref. [^, the valence quarks carry 47% of the total 
pion momentum at = 4GeV^. Note that, as one might expect. Fig. ^ also illustrates 
the fact that reproducing the momentum fraction is not sufficient to accurately determine 
the full shape of the distribution functions. From this point of view the agreement of 
the NJL evolved valence quark distribution with the SMRS92 parametrization is not 
entirely trivial. 



For comparison, let us also mention that early lattice calculations of Ref. [^, 
provided = 0.64 ± 0.10 scale ^ 4.84 ± 2.2GeV^ A recent and more accurate 
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Figure 6: Valence u-quark kaon/pion ratio of LO and NLO distribution functions in the 
NJL model at = 4GeV^ compared with phenomenological analysis. We take a total 
valence momentum fraction (xV")^ = {xV)k = 0.47 at = 4GeV^. Experimental data 



from Ref. 45 



lattice QCD calculation ||3 
0.56 ± 0.02 at the scale ^ 



extrapolated to the chiral limit yields the number (xV^) = 
5.8GeV^, a larger value than suggested by phenomenology 
3^, ^ and expected from a quenched approximation. The transverse lattice calculation 
of Ref. m gives = 0.86 ± 0.02 at ^ IGeV^ whereas that of Ref. ^ provides, 

still at very low scales ~ 0.4GeV^, a form for the distribution amplitude surprisingly 
close to the asymptotic value, 6x(l — x). From their parton distribution function one gets 
(xK) ~ 0.76. 



3.3 Kaon and Eta structure functions 



For the kaon and eta, we assume the same Qq as for the pion. For the K~^, this immediately 
leads to 



{x{uk-uk + sk-sk)) = 0A7 at = 4GeV^ . 



(9) 



Our LO and NLO evolved results for the parton distributions are shown in Fig. (Q). 
A practical parametrization of the corresponding initial condition may be found in the 
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Figure 7: Total valence Tr"*" and LO and NLO distribution functions in the NJL model 
at Q2 = 4GeVl We take {xV)^ = {xV)k = 0.47 at = AGeV\ For n+ we define 
V = u — u + d — d and for we have V = u — u + s — s. 



Appendix. Similar to the pion case, there are only small differences between LO and 
NLO evolution. The only known information regarding K structure functions is the ratio 
between the valence up quark distribution in the kaon and the pion, which was originally 
reported in Ref. MB and has been reanalyzed in Ref. p9|. In Fig. (0) we show the NJL 



results, together with the data obtained from Ref. Besides the LO result, already 

shown in our previous work [E^], we provide the NLO ratio, which does not differ much 
from the former and is in fair agreement with the experimental data. For the meson 
the momentum fraction for the up and strange valence quarks turn out to be 



(x 



{uk - uk)) = 0.20 



[X 



[SK - sk)) = 0.27 



at 



= 4GeV^ 



(10) 



As could be anticipated from Fig. (^), the difference for these momentum fractions between 
LO and NLO evolution are small and do not show up within the presented accuracy. 

Although a phenomenological analysis of the r] partonic distributions seems unlikely, 
for the sake of completeness we show in Fig. (|]) our results for the r] meson. We do this 
by evolving from the scale where a = 1.89(1.49) at LO (NLO) to = 4GeV^ the NJL 
distributions conveniently parametrized in the Appendix. As explained in our previous 
work E^, our description relies on a very particular ansatz which provides flavor mixing 
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without quark mass mixing. For the momentum fractions, we obtain 

{xu^) = (xdri) = 0.10 {xsri) = 0.08 at = AGeV^ . (11) 

As we have noted, the differences in parton distribution functions for massless and 
massive pions are tiny. In fact, even for K and r], many of the distributions are close to 
those of the massless pion. By comparing Fig. (|l|). Fig. (|) and Fig. (H) we point out 
the strong similarities in the gluon parton distributions between the 7i,K and r] mesons. 
Likewise, we also find very similar shapes for the sea distributions in the vr and K mesons, 
see Fig. (0) and Fig. (|^), as well as in the total valence distributions, see Fig. (|^. This 
is in agreement with having identical total valence momentum fractions for the pion and 
the kaon at = 4GeV^ 



4 Conclusions 

In the present work, we have computed the parton distribution functions of the lowest 
pseudoscalar mesons, namely vr, K and 1]. To this end we have used the Nambu-Jona- 
Lasinio distribution functions at the low resolution scale found in our previous work. In 
common with state of the art calculations, the corresponding sea and gluon distribution 
functions vanish at that scale, and are dynamically generated through standard GLAP 
evolution to higher Q^-values at LO and NLO approximations. For both vr and K we have 
assumed that the valence quarks carry 47% of the total momentum fraction at 4 GeV^. 
Despite the fact that a{Qo)/7[ ^ 0.5, the differences between LO and NLO evolution are 
small. The agreement between the u-quark valence distribution in the pion in the NJL 
model at LO and the phenomenological analyses is not spoiled at NLO. In addition, we 
have confirmed at NLO the successful description at LO of the ratio of the valence up 
quark content in the kaon with respect that of the pion. This provides one with some 
confidence in the validity of this approach to the study of structure functions, or GPD's 
in general. 

We have also presented LO and NLO sea and gluon distributions of the pseudoscalar 
mesons. For the pion, we find disagreement with the phenomenological expectations; the 
gluon and sea distributions come out to be softer in the high-x region and harder in the 
low-x region than the experimental analysis suggests. We have also provided results for 
the rj meson, which interest seems only theoretical, given the lack of experimental data. 
Our analysis, however, reveals some clear trends: all gluon distributions look strikingly 
similar, and the total valence tt and K distributions do not differ much. We hope these ob- 
servations to be useful to get further insight and guidance into the theoretical description 
of the poorly known meson structure functions. 

Appendix 



The 71, K and t] structure functions found in Ref. [|22[ may be conveniently written as 



u^{x) = 4(1 -x) = AN.gl^^-^ \mlFuuiml,x)] (12) 



''dm 



13 



uk{x) = sk{1 - x) 

Uri{x) = Uri{x) = drj{x) = drj^x) 



d 



m]^Fus{m\,x) 



-1 



q2 

\ JTJUU 



+ 



d 



iN, 



c 2 

9^ 



+ 



r]uu 



q2 



g2 



dm^ 

d 

dm'^ 



^T^'qFuui^f^rji x") 



(13) 
(14) 
(15) 



in the interval < x < 1. The Pauh-Villars regularized one-loop integrals are defined, 

1 



Fapip ,x) 



167r2 



J2 Ci log -x(l - x)p'^ + (1 - x)Ml + xMl + K^- 



(16) 



where Y.i Cif{Af) = /(O) - /(A^) + A^/'(A^). All other distribution functions are exactly 
zero, since we do not have gluons or sea quarks in the model. The meson-quark-quark 
couplings are defined in terms of the residues of the poles in the q — q scattering amplitude, 
and have the precise form needed to ensure the normalization conditions 

{u^{x)) = {d^ix)) = 1 
{uk{x)) = (sx(x)) = 1 

{Ur, + dn + Sr,) = 1 

{Ur, + drj + S^) = 1 . (17) 

The function Fq,^ satisfies the symmetry relation Fapip^ , x) = Fpa{p'^, l — x). This feature, 
along with the normalization condition, ensures the momentum sum rule. For the kaon, 
for example, one obtains 



{xuk{x) + xsk{x)) 
{xuk{x) + (1 — x)uk{x)) 



{xuk{x) + xuxi^ — x)) 
{uk{x)) = 1 . 



To apply the evolution method employed in Ref. |^ some analytical formula for 



the moments is needed. To obtain an approximate analytic formula for the moments, 
we note that for < x < 1, a convergent Taylor expansion of x dependence of Eq. (|T^) 
exists. Thus, the distribution functions may be accurately approximated by an nth degree 
polynomial, and the accuracy may be increased by keeping higher order terms. The pion 
and kaon distribution functions at Qq are accurately represented by 

{x, Ql) = dn+ (x, Ql) = 0.9535 + 0.2664a; - 0.2074a;2 

-0.1046a;^ + O.OlQOa;^ + 0.0400x^ - 0.0133a;^ (19) 
uk+{x,QI) = 1.1039 + 1.8071X- 1.0739x2 

-16.2227x3 + 33.1781x^ - 25.5372x^ + 7.1872x^ (20) 
SK+ix,Ql) = 0.4425 + 0.8593X -I- 1.7623x2 

-4.8611x3 + 13.2997x^ - 17.5858x^ + 7.1872x^ . (21) 



14 



In each case, the remaining quark and gluon distribution functions are assumed to be zero. 
For the r] meson, the expansion does not converge rapidly because one of the expansion 
parameters is 



M^x{l -x) ^ M2 



M2 



4M2 



0.8 . 



[22) 



The convergent series for Urj{x) we find to be given by 

'M2 + A2\ A2 



Ur,{x) 



A.,, 



In 



M2 



M2 + A2 



A2 



n=l 



(23) 



where a = M^/M^, /? = M^/{M^ + A^) and A„ 
for Srj{x) with the replacements M„ — ^ and 



0.09077. The same expression holds 
As = 0.36307. Although this series 
could be rearranged into a polynomial in x, it is easier to express the moments in terms of 
Euler complex Beta functions. In practice, 30 terms in the expansion are kept, providing 
a reasonable 0.1% accuracy. 
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